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A THEORETICAL AND EXPERIMENTAL STUDY OF THE NONLINEAR 


FLEXURAL VIBRATIONS OF THIN CIRCULAR RINGS 

By David A. Evensen 
Langley Research Center 


SUMMARY 


The nonlinear flexural vibrations of thin circular rings are analyzed by 
assuming two vibration modes and then applying Galerkin's procedure on the 
equations of motion. This procedure results in coupled nonlinear ordinary dif- 
ferential equations with time as the independent variable. The applied loading 
is taken to be harmonic in time, and approximate solutions to the equations are 
obtained by the method of averaging. One such solution involves the vibration 
of a single bending mode; a subsequent stability analysis shows that this 
single-mode response is valid only for certain combinations of amplitude and 
frequency. For example, when the driven mode exceeds a critical amplitude, non- 
linear coupling causes its companion mode to respond and participate in the 
motion. Approximate solutions are obtained for this coupled-mode case, and 
their stability is examined. The steady- state response curves contain an 
unusual "gap," where both the one- and two-mode solutions are unstable. These 
results were confirmed on an analog computer, and nonsteady vibrations were 
observed in the gap region. 

An experimental study of the problem was also conducted. Theory and 
experiment both indicate a nonlinearity of the softening type and the appear- 
ance of the companion mode. Measurements of the steady-state response are in 
good agreement with the calculated values, and the experimentally determined 
mode shapes agree with the form of the assumed deflection. 

The analytical and experimental results exhibit several features that are 
characteristic of nonlinear vibrations of axisymmetric systems in general and 
of circular cylindrical shells in particular. 


INTRODUCTION 


Current design of launch vehicles relies heavily on the use of thin-walled 
cylindrical shells as the primary structure. During powered flight, these 
cylindrical structures are often caused to vibrate to large amplitudes, in 
response to their environment. This problem has given rise to a number of 
theoretical studies of the nonlinear vibrations of thin cylindrical shells 



(refs. 1 to j). However, recent experimental results (ref. 4) suggest that 
such vibrations are still not fully understood. These studies prompted the 
investigation of a simpler but related problem, namely, the nonlinear flexural 
vibrations of a thin circular ring. 

Flexural vibrations of circular rings were originally analyzed in 1871 by 
Hoppe (ref. 5)j who presented the linear vibration frequencies and mode shapes. 
Shortly thereafter, Rayleigh (ref. 6) obtained similar results by using the 
approximation that the midsurface of the ring was inextensional. Since that 
time, the effects of midsurface extension, shear deformation, and rotary inertia 
have been investigated (refs. 7 to 9)> 

The first study of the elastic, nonlinear flexural vibrations of rings 
appears to be the 1959 work of Federhofer (ref. 10 ), who analyzed the free- 
vibration problem. The same problem was examined by Shkenyev (ref. ll) , who was 
primarily concerned with the dynamic stability of rings, as were Goodier and 
Mclvor (ref. 12). No experiments on the nonlinear vibration of rings have been 
reported in the literature, although linear vibrations have been studied exper- 
imentally (refs. 13 and 14 ) . 

The purpose of the present work is to study the forced, nonlinear flexural 
vibrations of thin circular rings. Only vibrations in the plane of the ring 
are considered, and the stress-strain law (Hooke's Law) is assumed to be linear. 
The nonlinearities examined here are geometric in nature and arise from the 
nonlinear terms of the strain-displacement relations. The ring is assumed to 
be of uniform rectangular cross section and to be relatively thin. With these 
restrictions, a detailed theoretical and experimental study of the problem was 
conducted. 

Agreement between theory and experiment was generally very good. The cal- 
culated response curves compare favorably with the experimental results, and 
measurements of the actual vibration form confirm the choice of the mode shapes 
used in the analysis. Both theory and experiment exhibit coupled vibrations 
that involve a driven mode and its companion mode. The results demonstrate 
several features that are characteristic of nonlinear vibrations of axisym- 
metric systems in general and of circular cylindrical shells in particular. 

This research was conducted while the author was a graduate student at the 
California Institute of Technology; the present report is a condensation of 
reference 15 • 


SYMBOLS 


The units used for the physical quantities defined in this paper are given 
both in U.S. Customary Units and in the International System of Units (SI) 

(ref. 16). Appendix A presents factors relating these two systems. 
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A 0 (t) 

A-n('t)^ B n (t ) 
A(t),B(t) 

Ag(t ) 

A 0 ,B 0 and 

b 

D 

d 

E 

e 

F 

F n 

G n 

h 

i* J 

K,Q 

k 

Mx>^xy>^y 
N x , ®bcy ) %■ 
n 
P 

<l(x,y,t) 


amplitude and time -varying part of uniform contraction, in. (cm) 
(see eqs. (8) and (9)) 

generalized coordinates associated with cos 22L and sin 2^ 
modes, in. (cm) (see eqs. (8) and (10)) 

slowly varying amplitudes, nondimensional 

C j. ( t ) expe riment al amplitude s 

$ 0 ,\|r 0 steady- state values of amplitudes and phases, respectively 
(see eqs. (F2)) 

ring length, in. (cm) 

Eh^ 

bending stiffness, -> r, in- lb (m-N) 

12 (l - v 2 ) 

displacement, in. (cm) (see fig. l4) 

Young’s modulus, lb/in^ (w/m^) 
base of Napierian logarithm 

load per unit length applied to ring in experiment, lb/in. (N/m) 
magnitude of assumed loading, lb/in. (N/m) (see eq. (l?)) 
dimensionless force 
ring thickness, in. (cm) 
integers 

Mathieu equation coefficients 

spring constant of drive wire, lb/in. (N/m) 

resultant moments per unit length, in-lb/in. (m-N/m) 

resultant forces per unit length, lb/in. (N/m) 

number of circumferential waves (n = 2, 3, 4 , . . . ) 
magnitude of load experimentally applied to ring, lb (N) 
radial load applied to cylinder 
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I ■■ I ■ ■■ 


R 

r 

t 

U 

u,v,w 


x,y,z 


Y,Z 
Pc, Ps 

A = ijr - 0 
&o 


e 


radial load applied to ring 
mean radius of ring, in. (cm) 
fraction of half wavelength 
time, sec 

stability variable used in equations (D7) 

displacement of a point on middle surface in axial, tangential, 
and radial direction, respectively 

coordinate in axial, circumferential, and radial direction, 
respectively 

Mathieu equation variables 

percent critical damping in cos ^ and sin ~ mode, 
respectively 

amplitude of shaker displacement, in. (cm) 
dimensionless nonlinearity parameter. 



e c 


corrected nonlinearity parameter (see eq. ( 29 )) 



e .€ direct strains 

e xy , e xz , e yz shear strains 

yiy 

£ c ( t ),£ s (t) nondimensional generalized coordinate associated with cos — 

ny 

and sin -r~- mode, respectively 


A,Ai 





v 


nondimensional stability variables 


Poisson's ratio 
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s c ( t )>s b ( 

p 

°xx.» °yy> 0 

CT xy 

r 

®(t ),lK T ) 

X = At + 
ft 

CD 

CD C 

“L 

(DM 

CD S 

V 2 

v* 

0(«2) 


t) perturbation in response of cos and sin mode, respec- 

R R 

tively (see eqs. (Dl) and (D2)) 
mass density, lb-sec 2 /in^ (kg/m3) 
zz normal stresses, lb/in^ (N/m^) 

shear stress, lb/ in^ (N/m2) 

nondimensional time, cD^t 
experimental phases 


slowly varying phases, nondimensional 


nondimensional frequency, cd/cd^ 
vibration frequency, radians/sec 

linear frequency of cos 2^; driven mode (see fig. 15 ) 

R 


nth linear vibration frequency, 
radians/sec 


E 

n 2 (n^ 

1 

H 

ro 

5 

~K> 1 

_pR 2 

(n 2 + l) 12(2 



1/2 


approximate vibration frequency, 
radians/ sec 


JL (n 2 - l) /h \ 2 

pR 2 12 (l - v 2 )' R '_ 


. l !/2 


linear frequency of sin 2 Z; companion mode (see fig. 15 ) 

R 


Laplacian operator, 


d 2 + 5 2 


etc 2 By 2 

biharmonic operator, V 2 ^ 2 

p 

order e 
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Subscripts: 
cr critical 

max maximum 

A dot over a quantity indicates differentiation with respect to time; a 
bar over a quantity indicates an average value. 


THEORY 


In this report, the nonlinear vibrations of a thin ring are analyzed by 
approximate techniques. The governing equations of motion are presented first, 
along with the associated boundary conditions. These equations are then reduced 
to ordinary differential equations by applying Galerkin's procedure. Two 
coupled nonlinear equations result, and approximate solutions to them are 
obtained by the method of averaging. The stability of these solutions is exam- 
ined, and improvements in the analysis are discussed. 


Equations of Motion 


The nonlinear 
(see appendix B) 


equations of motion of a thin ring can be written in the form 


^y 

Sy 



(la) 


D 



+ 




+ 






ph = q(y, t) 

at * 2 


(ib) 



Figure 1.- Ring geometry and coordinate system. 


where the circumferential force is given by 


% = Eh 


5v w l/5w\ 

Sy R 2 W/ 


(2) 


and the radial displacement w, the tangen- 
tial displacement v, and the applied load 
q are taken to be functions of only the 
circumferential coordinate y and the time 
t. (See fig. 1.) These results are 
obtained by a specialization of the analo- 
gous equations for thin cylindrical shells. 


Boundary conditions . - Since the ring 
is circular, the displacements must satisfy 
the following continuity requirements: 
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w(y,t) = w(y + 2*R,t) 

^y,t) - + atR,t) 

oy dy 1 - 


(i = 1, 2) 


(3a) 


v(y,t) = v(y + 2rtR,t) 


bv 

dy 


(y,t) =■ |^(y + 2jtR,t)j 


(3b) 


in addition to the equations of motion. From a physical standpoint, these 
boundary conditions insure that the displacements, slope, and other related 
variables are continuous in the circumferential direction. 


Simplified equations of motion .- Previous work (ref. 15 ) has shown that for 
flexural vibrations, equations (l) can be simplified considerably without losing 
the essential features of the problem. This simplification is accomplished by 

(1) assuming that the midsurface circumferential strain is zero 

(2) neglecting the effect of tangential inertia 
The first of these conditions requires 


I _ dv , w , l/dw\ _ 0 
z=0 by R 2ydy J 


( 4 ) 


where the nonlinear strain-displacement relations (appendix B) have been used. 
Similarly, neglecting tangential inertia in equation (la) results in 



dy 


0 


(5) 


Combining equations (2), (^), and ( 5 ) gives Ny = 0. Using this result in 
equation (lb) yields the simplified equation of motion for w: 


D 



+ 





+ ph 



q.(y>t) 


(6) 


Note that the continuity conditions on w and v are not altered by the 
preceding discussion. 


With these simplifications, the remaining problem is to solve equation (6) 
subject to the continuity and inextensionality constraints. Since the latter 
requirement (eq. (4)) is nonlinear and obviates exact solutions to equation (6), 
it is necessary to resort to approximate techniques. 
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Nonlinear Ordinary Differential Equations 
for Vibration of a Thin Circular Ring 

Approximate solutions to equation (6) can be obtained by assuming the shape 
of the deflection in space. This approach is commonly used in nonlinear vibra- 
tions of structures (cf. refs. 17 to 19) and reduces the problem to one 
involving nonlinear ordinary differential equations in t . 


Choice of deflection modes .- The most general radial deflection compatible 
with the continuity requirements is 


v(y,t) 


n=0 


A n (t)cos ^ + B n (t)sin 
R 


ny 

R 


(7) 


where A n (t) and B n (t) are periodic in time. This case has been analyzed in 
detail in reference 15, but the majority of the results can be obtained by using 
the following two-mode approximation: 


v(y,t) = A n ( t ) cos SC + B n (t )sin SC + A 0 (t) (8) 

ny ny 

Here cos — and sin — are the linear vibration modes of the ring; and since 
R R 

only flexural motions are considered, equation (8) is restricted to cases where 
n 2. The n = 0 mode is (by itself) an axisymmetric motion involving 
stretching of the midsurface, and the n = 1 modes correspond to displacement 
of the ring as a rigid body. The Aq term in equation (8) is related to A n 
and B n by the inextensionality constraint. 

Determination of A 0 (t) in terms of A n (t) and B n (t).- Solving equa- 
tion (lj-) for dv/dy and substituting for w from equation (8) gives 

8v _ _ w l/dw\ 

8y R 2 y8y J 


Ao("t ) 
' R 


p O 

A n (t) + B n (t) + Terms periodic in 


y 


To satisfy the constraint v(y,t) = v(y + 2?cR,t), the terms that depend solely 
on time must be equated to zero. This procedure yields 


A 0 (t) = - 




(9) 
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and equation (8) can then he rewritten as 


v(y,t) = A n (t)cos + B n (t)sin ~ 



( 10 ) 


for the assumed radial displacement. This expression is compatible with the 
continuity and inextensionality constraints and can be used with Galerkin’s 
procedure to satisfy approximately the simplified equation of motion for w. 


App lication of the Galerkin method .- To apply Galerkin’s procedure, equa- 
tion (10) is substituted for w in equation (6). The resulting expression is 
then multiplied by a weighting function associated with A n (t) and integrated 
with respect to y from 0 to 2 jxR. This procedure yields an ordinary differ- 
ential equation involving primarily A n (t). An equation for B n (t) is obtained 


in a similar fashion; both equations are coupled in the nonlinear terms, 
weighting functions used are 


and 


8w 

ny 

n 2 A n 

8A n 

= cos — “ 

2R 

8w 

= sin 2 - 

n2B n 

SB n 

R 

2R 


The 


(11a) 


(lib) 


respectively. 

Carrying out the operations just indicated and then nondimensionalizing 
the results yields the following coupled equations: 


d ^c 

dT 2 


d ^c „ e ^c 

+ Vc + t c + ~f 


dT 


2f / j ► \ 2 




S &T 2 \ ^ / 


= G n cos fir 


d 2 ? 


s J. OQ + £ + fii 

j c: 


dT 


2 + 2Ps dT " s 


s . + («-f + c 

S dT 2 V dT ) 


d ^C , 

C dT 2 \ dT 


= 0 


(12a) 


(12b) 


where the dimensionless variables are 


5 = ^ 

s c h 


t = ^2; 
h 


= 


n = JH. 
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e 



G n = 


P 


xi 


jtRph 2 ^ 2 


and 

(l) for simplicity, the loadixig has been taken to be 

F 

q(y,t) = -§■ COS ^ COS OJt 
JtK x\ 


(13) 


(2) two minor nonlinear terms have been discarded 

(3) modal damping terms have been inserted to study the first-order effects 
of small viscous damping 

Although these results were derived by a Galerkin procedure, they can also 
be obtained by the Rayleigh-Rit z method. In general, these two approximate 
techniques are not equivalent; they can be made to coincide, however, by the 
proper choice of weighting functions (ref. 20). It was for this reason that 
the weighting functions and c3w/c)B n were used to obtain equations (12). 


Approximate Solutions to Nonlinear Equations 


The previous section indicated the manner in which the equations of motion 
are reduced to ordinary nonlinear differential equations. These equations still 
cannot be solved exactly, but approximate solutions to them can be obtained by 
the "method of averaging." (For a discussion of this technique, see ref. 21, 
ch. V.) 


Vibrations involving a single bending mode ^cos are considered first, 

and nonlinear single-mode response curves are presented. When the stability of 
the single-mode solution is examined, the vibrations are found to be stable for 
only certain combinations of amplitude and frequency; this result makes it 

ny , ny 

necessary to examine vibrations where both cos — and sin — participate m 

K R 

the motion. The method of averaging is used to obtain approximate solutions 
for this two-mode case and also to examine their stability. 


Response of a si ngle bending mode .- Inspection of equations (12) reveals 
that possible solutions are £ s (t) = 0 and £ (t) j= 0, where the latter 

satisfies 


a 2 e. 


dr" 


+ 2(3, 


dT 


+ 




d 2 £ 

t C 

3 r» 


dT' 


+ 

2 \ dT 




= G n cos Gt 


( 14 ) 
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To obtain approximate solutions to equation (l4) by the method of averaging, let 

£ c (t) = A(t)cos[at + <D(t)| (15) 

where A and 3> are presumed to be slowly varying functions of t. When this 
expression for £ c is substituted into equation (l4) and the appropriate aver- 
ages are carried out, two equations result (see appendix C). 

\ 2—3 

(l - n 2 )A - = G n cos ¥ ( 16 a) 


-20 C £2A = G n sin $ (l6b) 

Here A and $ are average values (over one period) of A(t) and $( t ) . 
Squaring and adding equations (l6) results in one equation involving X: 


2 

+ 4p c Va 2 = G n 2 (17) 

For given values of G n , e, 0 C , and ft, equation ( 17 ) can be used to compute 

A. Then equation (l6b) can be solved for $, and the approximate solution to 
equations (12) becomes 


(l - d, 2 )a - 


2-3 

eO 


C c (t) = A cos(Gt + $) 


( 18 a) 


£ s (t) = 0 (18b) 

for vibrations where only one bending mode ^cos responds. A typical reso- 

nance curve is shown in figure 2, which demonstrates the nonlinearity of the 
softening type. 

The case of free vibrations may be obtained by putting G n and 0 C equal 
to zero in the preceding results. Equations (l6) then yield 

fl = 1 - + o(e 2 ) ( 19 ) 

8 

which is the so-called "backbone curve" for free nonlinear vibrations. This 
curve is illustrated in figure 3 tor various values of e . 

Stability of the one-mode response .- The stability of the preceding solu- 
tion was investigated by perturbing £ c (t) and £ s (t). A study of the 
resulting Mathieu equations indicates that within order e2 
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Figure 2.- Typical single-mode response. Dr 
e = 4.2 x 10" 4 ; p c = 2 x io _3 ; G 

(1) perturbations of £ c 

1 

( 2 ) perturbations of £ s are unstable within the region 

eA 2 eA 2 

o o 


iven mode, cos 
n = °- 10 - 


Figure 3.- Backbone curves for various values of £. 
Nonlinearity parameter, £ = 


are unstable within the area bounded by 

-2 


-2 

3eA 

8 


< n < l 


eA 

8 


( 3 ) both types of perturbations are also unstable 


2 = 


111 
2 " 3 ' h’ 


in narrow regions near 


Details of the stability calculations are given in appendix D, along with 
an interpretation of the stability boundaries. The results of the stability 
analysis (for the vibration of a single mode) are summarized in figure 4. In 
one region, the solution is unstable with respect to perturbations of and 


12 



jumps occur in the response. In an 
adjacent area, perturbations of £ s are 

unstable, and the companion mode begins 
to vibrate; this case of coupled vibra- 
tions is analyzed in the following 
section. 


Re spon se of the self- coupled bending 
modes .- When £ b (t) and £ c (r) both 

oscillate, their coupled vibrations are 
governed by equations (12). As in the 
one-mode case, the method of averaging is 
employed to obtain approximate solutions. 
To use this method, let 


£ c (t) = A(t)cos[7h + *( T g 
£ s (t) = B( T )sin j^h T + *( T j] 


( 20 ) 


where A, B, 0, and r|r are presumed to 
to be slowly varying functions of t. 


Substituting the expressions for 
and £ s from equations (20) into equa- 
tions (12) and applying the method of 
averaging gives 



Figure 4.- Stability diagram for single-mode response. 



20A 


d$ 

dT 


eA 

4 


+ | nA 2 M 

2 dT 


P — 2 — 

il^B cos 2A 


fiB dB 
2 dr 


sin 2A 


+ 2B 2 fl - ~ cos 2 a]|^ 


= G n cos 0 


(21a) 




.20 fi - 2BM - TVS 


sin 2A - ~ cos 24 §2. 

2 dT 


+ sin 24 I = G n sin 0 
ar / 


(21b) 


(l - A 2 )b - 2CIB ^ -fl^cos 2A + VpB 2 + VA 2 (l - J cos 24^ 

dT 4 \ 2 /d T 


+ ££ sin 24 | 2B 2 

2 dT 2 dT 


= 0 


(21c) 
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2G S® + 23 S 2B - —-L 2 ! 2 sin 2A - — — 
dr s 4 \ 2 d T 

+ — cos 2A — + sin 2A — ) = 0 (21d) 

2 dT 2 dr / 

where A(t), B(t), $(t), and i|/(t) have^been replaced by their average 
values over one cycle and are denoted by A , B, i, and if. The average 
phase difference A is 

A = t - ® (22) 

For steady-state vibrations, the average values _A, B, <±>, and i|r are 
all constant; in this case, their time derivatives (dA/dT, dB/dT, etc.) are all 
zero, and equations (21) reduce to 


(l - G 2 )a - £^(a 2 - B 2 cos 21 ) - G n cos $ 

(23a) 

-23 C GA - A B 2 sin 2A - G n sin $ 

( 23 b) 

(l - G 2 )b - ^^(b 2 - A 2 cos 2 a) - 0 

(23c) 

23c.GB - A 2 sin 2A - 0 

S 4 

(23d) 

Equations (23) may be solved simultaneously for A, B, 3>, and 
approximate solution to equations (12) then becomes 

i|r, and the 

£ c (t) = A cos(Gt + i>) 

(24a) 

£ b (t) = B sin(GT + t) 

(24b) 


Details of the_calculations are_given in appendix E. Response curves of the 
variation of A with 0, and B with 2 are shown in figure 5, where the 
particular values 


G n = 0.1 (Force = G n cos Gt) 

3 C = 3 S = 2 x 10“^ (Viscous damping) 


l4 



and 


e 



= 4.2 x io -1+ 


(Nonlinearity parameter) 


were used. 

A surprising feature of the two-mode response was the appearance of a 
"gap," or discontinuity in the solution. Sl.ightly to the left of the resonance 
peak in figure 5, the approximate solutions ifor £ c and £ s break down. The 
results of a subsequent stability analysis showed that the gap in the response 
coincides with a narrow region in which both the one-mode solution (eqs. (l8)) 
and the two-mode solution (eqs. (24)) are unstable. The experiments suggest 


that a "beating" response exists in 
this area, with the ring vibrating 
first in one mode, then in two modes, 
then back to one, and so forth. 

Analog- computer studies (dis- 
cussed in a subsequent section) veri- 
fied the calculated response curves, 
and nonsteady vibrations were found to 
occur in the gap region. Similar 
responses involving gaps have been 
observed in the case of fuel sloshing 
and in other nonlinear problems. (See 
refs . 22 and 23 . ) 

Another curious result demon- 
strated in figure 5 is that in some 

cases the companion mode (sin which 

V R \ 

is not driven by the forcing function ] 
can vibrate to larger amplitudes than 
the driven mode (cos Responses of 

this type occur for ft < 1 and were 
detected experimentally. 

It will also be noted that in the 
two-mode case the amplitude of the 
driven mode is considerably reduced 

from what it would be if sin — did 

R 

not vibrate. A comparison of the 
driven-mode response for both cases 
(one mode and two modes) is shown in 
figure 6. 



(a) Companion mode, sin 



(b) Driven mode, cos-^-. 

Figure 5.- Typical coupled-mode response, e = 4.2 x lCf^; 
Pc = P s = 2 * 10‘ 3 ; G n = 0.10. 
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Finally, both figures 5 and 6 show that the two-mode response is a lmost 
linear, at least in those areas where the vibrations are stable. 


Stability of the coupled-mode response .- As for the one-mode case, the sta- 
bility of the coupled-mode vibrations can be found by perturbing the steady- 
state response. _The resulting analysis shows that within order real non- 

zero values of B do not exist for 

Q > 1 + 

8 

and the two-mode solution is unstable for 


a < i - 

8 

These results are presented in figure J, and the corresponding stability cal- 
culations are outlined in append!:*: F. 

A comparison of the preceding results with the results for the unstable 
regions given previously shows that the one-mode and the two-mode solutions are 
both unstable in a narrow region, given by 



two modes, e = 4,2 * 10‘ 4 ; = p s = 2 x i(r 3 ; G n = 0.10. 



Figure 7.- Stability diagram for coupled-mode response. 


16 



1 


o o 

As noted previously, this region coincides with the gap in the two-mode 
response. The gap in the solution is apparently caused when the assumption of 
slowly varying amplitudes and phases inherent in the theory no longer applies 
to the actual motions of the ring. 

This discussion concludes the one- and two-mode analysis in this report. 

The calculations can be extended to account for the effect of tangential inertia 
and additional nonlinear terms which have been previously neglected. The gen- 
eral case, involving an infinite number of bending modes, has also been ana- 
lyzed. (See ref. 15*) These investigations indicate that the basic features 
of the vibrations are adequately given by the one- and two-mode study just pre- 
sented. The results of the more refined approaches are briefly outlined in the 
following section; details of the calculations may be found in reference 15. 


Improvements in the Analysis 

Throughout the present study, many simplifying approximations have been 
used. In order to determine when these approximations are valid and when they 
break down, it is necessary to consider some refinements in the analysis. These 
refinements yield significant improvements in the theory and bring it in closer 
agreement with the experiments. 

Discussion of tangenti al inertia and other effects .- Since the effect of 
tangential inertia has been neglected, equation (6) yields a slightly inaccu- 
rate result for the linear vibration frequency. Including the tangential 
inertia term (eq. (la)) corrects this discrepancy, and repeating the analysis 
gives results in the form of equations (l 6) and (19), where ft is replaced by 


cd Frequency of vibration 

<*>L Linear vibration frequency 


and e is replaced by 


e i 



(25) 


In other words, the form of the results is unchanged, but the parameters ft and 
e are altered to account for the tangential inertia. 


Other small corrections, such as those due to shear deformation, are gen- 
erally unimportant if the ring is sufficiently thin. For instance, Buckens 
(ref. 8) has shown that shear deformation and rotary inertia may be neglected if 


n 2 h 2 


3R 


« l 
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This result Is in agreement with the analogous result for beams - namely, that 

p 

shear effects can be neglected if the (depth/length) ratio is much less than 
unity. 

Similarly, the effect of extension of the midsurface of the ring is usually 
unimportant for bending vibrations oi’ thin rings. This result is well estab- 
lished for linear vibrations (see rei’s. 6, 8, and 12) and carries over in the 
nonlinear case as well (see ref. 15 in this regard). 


C orrections fo r additional nonli. nearities . - The preceding corrections are 
all relatively minor, and their influence is fairly well predicted on the basis 
of linear vibrations. A major improvement in the nonlinear analysis can be 
obtained when additional nonlinear te:rms are included in the strain-displacement 
relations. Up to this point, the calculations have all employed the usual non- 
linear expression for the midplane circumferential strain, namely. 


e 


yy 


z=0 


/ \ 2 

5 v + w + i/8w ) 

&.y R 2 \8yy 


(26a) 


which is commonly used in nonlinear shell studies (e.g., refs. 1 to 4 ). How- 
ever, equation (26a) is an approximation to the more complete strain-displacement 
relation (ref. 24 ) 


dv w 1 

e yy z=0 8y R 2 


8w 

Sy 


-iO + (i + D‘ 


(26b) 


When the deflection w(y, t) is assumed in the form (eq. (8)) 


w(y,t) = A n (t)cos ~ + B n (t)sin + A 0 (t) 


ny 


and equation (26a) is used in the inextensionality condition, A Q (t) is found 
to be (eq. (9)) 


Ao(t) = - S 


4 h 


A n 2 (t) + B n 2 (t) 


Then the deflection is (eq. ( 10 )) 


w(y,t) = A n (t)cos + B n (t ) sin ^ - 2 _^ n 2 (t) + B n 2 (t)_ 

and the coefficient n 2 / 4 R eventually becomes the nonlinearity parameter 

e = • Thus, the coefficient of |A n 2 (t) + B n 2 (t)j in the preceding 

expression for w influences all the previous results (e.g., the response 
curves, the backbone curve, and stability boundaries) by means of e. 
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In a similar fashion, when the more exact equation (26b) is used in the 
inextensionality condition, a more accurate expression for A 0 (t) results 
(see ref. 15 ) 

A-oU) = “ ~ -g-'j [A n 2 ( t ) + B n 2 ( t )J + Terms of order (2j) 

If the higher order terms are neglected, the improved version of the deflection 
then becomes 

v(y,t) = A n (t )cos ZjL + B n (t )sin SL _ j£^L _ J^ 2 [ An 2 (t) + B n 2 (t)] (28) 


Since this expression for w has exac tly the same form as the one which was 
previously used (eq. (lO)), the analytical results are altered only through the 

new coefficient of + B n 2 (^)J • That is, the results of the previous 

calculations can be modified to account for the additional nonlinearities by 

replacing e = 




This correction can be combined with the one which accounts for tangential- 
inertia effects by first defining a corrected nonlinearity parameter e c 


e c = 



(29) 


Then all the previous analytical results (response curves, etc.) may be cor- 
rected for both the additional nonlinearities and tangential inertia by 
replacing Q by cd/cd^ and e by the expression for e c (eq. ( 29 ))- The 

corrections are most noticeable for low values of n, with the effect of the 
additional nonlinearities predominating. For example, when n = 2, e c is 
about one-fourth as large as e (i.e., the nonlinearity is reduced by nearly 
a factor of 4). Values of e c /e for various values of n are given in the 
following table: 


Mode number, n 

e c/ € 

2 

0.261 

3 

• 571 

4 

• 734 

5 

.820 

10 

• 950 
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The bulk of the experimental results were obtained for the n = 4- mode, and 
the comparison of theory and experiment was made by using e c = 0. 73^-e • The 
experimental work is discussed in the following section. 


EXPERIMENTS 


Apparatus 

A thin seamless copper ring was used in the experiments; it had a radius 
of 4 inches (l0.l6 cm), a thickness of 5-1 X 10 - 3 inch (0.013 cm), and a 
length in the axial direction of 0.988 inch (2.51 cm). The ring was supported 
by four very thin suspension threads, equally spaced around the circumference, 
as shown in figure 8. 


Radial motions of the ring were measured by two inductance pickups which 
operated in a push-pull fashion. These deflection sensors were mounted on a 
fixture with a large bearing that allowed them to move circumferentially around 
the ring. Vibrations of the ring were excited by means of an electrodynamic 
shaker. The shaker was connected to the ring by a fine tungsten drive wire 
(0. 001- inch (0.003-cm) diameter), which served as a very soft coupling spring 
between the shaker and the ring. (A simplified sketch of the apparatus is shown 
in fig. 8.) This arrangement made it possible to estimate the amplitude of the 
force which was experimentally applied to vibrate the ring. 



Figure 8.- Simplified schematic of experimental apparatus. 


The response of the 
ring was analyzed, and the 
vibration modes were iden- 
tified by means of Lissajous 
figures . These and other 
technical details of the 
experiments are discussed 
in appendix G. 


Measurement of Mode Shapes 


The mode shapes were 
measured by exciting one 


mode cos 


ny 

R ’ 


n = 4) 


and 


recording the amplitude of 
the response at intervals 
along a half wavelength. 
Measurements were made sep- 
arately on two different 
half waves; the amplitudes 
ranged from about 1 to 27 
times the thickness of the 
ring. These amplitude 
measurements are shown in 
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figure 9> where the solid lines 

correspond to the deflection shape Assumed shape, sinwr 

assumed in the analysis . As indi- ° Ex P erimental data 

cated in figure S>, the main vibra- 

tion shape was virtually independ- - 0 " -0 _ 

ent of amplitude. Additional A =27ih 

results presented in reference 15 0 , , max , , 

include a plot of the root -mean- 
square response around the circum- 
ference of the ring and measurement 1 d^' 0 "” -0 0 

of the motions which occurred at 

the nodes of cos ^ for the 0 1 , max , 

^ A/A 

single-mode response. max 

Response Curves for the A max = 9 - 6h 

1 1 1 1 

n = 4- Mode 

Some typical experimental — 3 ° 

results are shown in nondimensional A = i.36h 

form in figure 10 for four differ- <U^ 1 , max , , 

ent values of the input force. The 0 2 A - 6 - 8 1.0 

magnitude of the input force was Fraction of half wavelength, r 

held constant for the individual 

response curves by maintaining a Figure 9.- Comparison of assumed mode shape and experimental results 

fixed displacement of the shaker for various amplitudes. 

(25 q , peak-to-peak) during each 

run. An electronic counter was used to measure the period of the forcing func- 
tion at each data point; by this procedure it was possible to determine the fre- 
quency ratio co/cdl very accurately. The amplitude of A n (t) was measured at 

one antinode of cos S2L, and the phase of the response 3>. was determined there 

with a standard phase meter; the magnitude of B n (t) was measured at an adjacent 

node of cos To facilitate comparison with the analysis, these experimental 

amplitudes were nondimensional! zed on the ring thickness. 


A 

= 27. lh 

max 


A 

- 15. 8h 

max 

1 

1 

A 

= 9.6h 

max 

1 

1 


A = I.36h 
max 


Fraction of half wavelength, r 


Figure 9.- Comparison of assumed mode shape and experimental results 
for various amplitudes. 


Corresponding theoretical results are shown by the solid lines in fig- 
ure 10. The backbone curves were calculated by using equation ( 19 ) with 

in place of and e c substituted for e; both charges were made to account 
for the influence of tangential inertia and the additional nonlinearities dis- 
cussed previously. The same corrections were used in computing the forced 
response, which was obtained from equations (l6) for the undamped case. The 
input force used for these calculations was computed from the shaker displace- 
ment and the spring constant of the drive wire. A detailed comparison of the 
calculated and measured response for the driven mode is given in figure 11, 
which is for a peak-to-peak shaker displacement of 28 Q = 0.3 inch (0.j6 cm). 
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Amplitude ... Amplitude 

Thickness' Thickness 



Figure 10.- Response curves for n = 4 mode. 


Amplitude . Amplitude 

Phase, $a Thickness’ Thickness’ 



(c) Shaker displacement: 

28 0 = 0.300 in. (0.762 cm). 

Figure 10.- Concluded. 


(d) Shaker displacement: 

25 0 = 0.400 in. (1.016 cm). 
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Discussion of the Experimental Results 

As shown in figure the experimentally determined mode shapes were 
nearly perfect sine waves along any half wavelength. This result was found to 
be independent of the maximum amplitude of the motion, as figure 9 indicates. 

For the single-mode response, the displacements at the nodes of cos 2^ were 

R 

found to increase linearly with the square of the maximum amplitude. (See 
ref. 15*) This behavior is in agreement with the analytical results. 

The response curves of figure 10 
show good correlation with the theo- 
retical backbone curves, and figure 11 
shows a frequency difference of less 
than 1 percent between theory and 
experiment. This slight discrepancy 
was probably due to small errors in 
estimating the force applied to the 
ring. 

The experimental response curves 
exhibited the jump phenomena and the 
appearance of a secondary resonance 
peak for the driven mode. The latter 
peak resulted from the fact that, in 
the experiment, the driven mode and 
the companion mode had slightly dif- 
ferent natural frequencies. Because 
of this disparity, it was not possible 
to make a quantitative comparison of 
the coupled-mode response with the 
corresponding theory. However, quali- 
tative agreement with the analysis was 
obtained for the two-mode case and 
included the result that in some 
instances the companion mode vibrated 
to larger amplitudes than the driven 
mode. One major feature of the theory which was not detected experimentally was 
the so-called "gap" in the response; the disappearance of the gap was apparently 
due to the difference in the natural frequencies of the coupled modes. The only 
results which might have been related to the gap were occasional nonsteady 
beating responses. These were sometimes observed to the left of the companion- 
mode resonance peak. 

In order to understand more fully the response of the ring in the gap 
region and to verify that such a gap did theoretically exist, equations (12) 
were solved on an analog computer. The analog solutions confirmed the previous 
analytical results, as shown in figures 2 and 5- Nonsteady vibrations with 
beats occurred in the gap region. The analog procedure is outlined in 
appendix H. 



Figure 11.- Calculated and measured response. n = 4 mode; 
26 0 = 0.300 in. (0.762 cm); e c = 3.08 x 10 -4 . 
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CONCLUDING REMARKS 


The nonlinear flexural vibrations of thin rings were analyzed by choosing 
vibration modes and applying Galerkin’s procedure. The vibrations were assumed 
to involve no stretching of the midsurface of the ring; this assumption was 
found to be adequate for the study of flexural vibrations. In the analysis, 
only one mode was directly driven by the forcing function; nevertheless, it was 
necessary to include two vibration modes in the calculations because under cer- 
tain conditions, nonlinear coupling caused the companion mode to respond and 
participate in the motion. In other cases, the single-mode response was 
sufficient . 

Significant improvement in the theory was obtained by including the effect 
of additional nonlinearities in the strain-displacement relations. Retaining 
the effect of tangential inertia also improved the calculations, but to a lesser 
extent. Both these modifications lose their importance as the mode number 
increases; however, they combined to decrease the nonlinearity parameter by more 
than 25 percent for the third flexural mode (n = k) . 

The experimental results were found to be in good agreement with the anal- 
ysis, both qualitatively and quantitatively. Theory and experiment both exhib- 
ited the jump phenomenon, nonlinearity of the softening type, and the appearance 
of the companion mode. The measured responses were in good agreement with the 
calculated values, and the experimental mode shapes demonstrated the appropriate- 
ness of the deflection form employed in the analysis. 

The results of the present study are characteristic of the nonlinear forced 
vibrations of axisymmetric elastic bodies. In such structures, the nonlinear 
forced vibration of one mode often results in the response of both the driven 
mode and its companion. Vibrations of this type occur because of the nonlinear 
coupling that exists between the modes involved. The nonlinear forced vibra- 
tion of thin circular cylinders, thin circular cones, and other thin axisym- 
metric structures can be expected to exhibit similar behavior. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August K, 1965* 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (Si) was adopted by the Eleventh General 
Conference on Weights and Measures, Paris, October i960, in Resolution No. 12 
(ref. 16 ). Conversion factors for the units used herein are given in the fol- 
lowing table: 



U.S. 

Conversion 


Physical quantity 

Customary 

factor 

SI unit 


Unit 

(*) 


Force 

lb 

4.45 

newtons ( N ) 

Frequency . 

cps 

1 

hertz (Hz) 

Length 

in. 

0.0254 

meter (m) 

Load 

Ib/in 2 , 

6895 

newt on/met er2 (ll/m^) 

Mass density 

lb-sec2/in b 

27 679.9 

kilogram/meter 3 (kg/m 3 ) 

Moment per unit length . . 

in- lb/in. 

4.45 

meter-nevton/meter (m-N/m) 

Stiffness 

lb/in. 

175.2 

newton/meter (N/m) 

Young ' s modulus ..... 

lb / in^ 

6895 

newt on/met er 2 (N/m 2 ) 


ji, 

Multiply value given in U.S. Customary Unit by conversion factor to obtain 
equivalent value in SI unit. 


Prefix to indicate multiple of units is as follows: 


Prefix 

Multiple 

centi 

10“ 2 
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APPENDIX B 


THE EQUATIONS OF MOTION 


The equations of motion for a thin circular ring become nonlinear when 
finite displacements are considered. The purpose of this appendix is to pre- 
sent a brief development of these equations. 


Useful Results From Shell Theory 

The equations for rings can be obtained in a straightforward manner by a 
specialization of the analogous equations for thin cylinders. The latter equa- 
tions have long been employed in nonlinear analyses of thin cylindrical shells; 
they are often derived (e.g., refs. 1 and 2) and are simply repeated here: 


bx 


3n. 


.2 

= P h ^ 


dN. 


xy 


by 

dN. 


c)x 


+ 2 ^xy + 3%- + b 

bx 2 bx by ^y 2 Sx 




by 

dw 

bx 


- = ph 


afv 


+ ir 


by 


N $1 + N 
% ^ + N xy 5x 


(Bl) 

(B2) 


N- 


R 


f + q.( x ^y^ t ) = pk — 


■W 


St" 


(B3) 


The variables u, v, and w are the midplane displacements in the axial, cir- 
cumferential, and radial directions, respectively. (See fig. 1.) Forces and 
moments per unit length (N x , Ny, Mjj, etc.), defined in terms of the stresses, 
are 


N x = 


M x = 


nh/2 

r h/ 2 


ph/ 2 

/ . *xx dz 

J -h/2 

N X y - / a^y dz 

J -h/2 

Ny = 

J- h/2 ^ ^ 

oh/2 

r h /2 


oh/2 

/ cr-^z dz 

; -h/2 

^xy - / CT xy z dz 

J -h/2 

% = 

/ Oyv z dz 

J -h/2 " J 


> (B4) 


The applied load q(x,y,t) acts in the radial direction. 

Equations (Bl) to (B3) contain both forces and moments as well as dis- 
placements. Before specializing these equations to apply to rings, it is 
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desirable to express them in a more workable form. To achieve the desired form, 
the following procedure is used: 

(1) The strains are related to the displacements u, v, and ' w 

(2) Hooke's law is used to relate the stresses to the strains 

( 3 ) By combining the relations from (l) and (2), the stresses are written 
in terms of the displacements and then employed in equations (Bl) to (B4). 

The strain-displacement relations are approximated by 


e xx 


8u l/dwY^ 8^w 

dx 2 \ox / g x 2 


"A 


e yy 


8v w + l/bv \ 

by R 2\8y J 



(B5) 


= ^ 2z — v 

6x ^ r 8x by bx by bx by^ 

Transverse shear deformation is neglected, and the transverse normal stress is 
assumed to be negligible. These assumptions may be written as 

e xz = 0 e yz = 0 a zz, = 0 

These relations are commonly employed in nonlinear analyses of thin cylinders. 

give rise to nonlinearities In the equations of motion. 


Terms like —I t— 


bw 

by 


Next, Hooke's law is used to relate the stresses to the strains and with 
the preceding assumptions yields 

E , 

°yy 1 _ v 2 \ e yy 

E / 

°xx pV £xx + 

1 - 

Gx y = 2(1 t v ) €x ^ 

where E is Young's modulus and v is Poisson's ratio. Employing these 
stresses and equations (B 5 ) in equations (Bk) and carrying out the integration 
through the thickness gives the following forces and moments in terms of the 
displacements: 


ve xx) 

V€ yy) 
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chi l/dw 
Bx 2\^x 


Hx= i^ + KtO + * 


bv + w + l/dw Y 
dy R 2\dy/ 


„ Eh (Sv 5u , 5w Sw 

N„„ = — -It— + — + 


x y 2(1 + v)\dx dy dx dy 


/.2 2 ' 

My = -D ^ + V ^ 

7 \dy 1 2 3 dx 2 / 


A * 2 .2 > 

ha -nl ° W O W 

M x = " D i— + v — 


(dx 2 dy £ 


Mxy - -D(l - v) 


w 


dx dy 


where the bending stiffness D = 


Eh" 


12(1 - v 2 ) 


(B6) 


Derivation of Ring Equations 

To specialize the cylinder equations to a ring, the following assumptions 
are made : 

(1) The displacements w and v, as well as the radial load q, are 
taken to be functions of only the circumferential coordinate y and time t. 

(2) The thickness and the width of the ring are both taken to be constant. 

The ring is assumed to be thin, so that ( — ] is negligible in comparison 
with unity . ' ' 

( 3 ) The forces N x and N X y are assumed to be zero throughout the ring. 
The boundary conditions require that N x and N X y vanish at the ends of the 

ring (at x = 0, x = b (fig. l))j for flexural vibrations, rapid changes in 
the x-direction are not anticipated, and it becomes reasonable to assume that 
N x = 0 and Nxy = 0 throughout the ring. 
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When these assumptions are used, equations (b 6) reduce to 


Ny — Eh 


N x = 0 


N xy 0 


*v-- D S 


Bv + w + l/Bw V 
By R 2 W/ 


2^ 


M x = -vD 


Mxy = 0 


By 

B 
By 2 


a 2 

o w 


(B7) 




Substituting equations (B7) in equations (Bl) to (B3) yields the following 
ring equations: 


where 


0 = ph 



SNy 

By 


ph 



cAw 

By 4 




+ ph ~~ = q(y,t) 
Bt 2 


N 


y 




+ 



(b8) 

(B9) 

(BIO) 


(Bll) 


Equation (B8) is unimportant for the problems dealt with herein and will hence- 
forth be dropped. 

Note that there is a similarity between equations (B9) and (BIO) and the 
following analogous equations for a vibrating beam: 


BT _ ~ B 2 v 


(B12) 
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El 




+ pS 


a 2 w 

St 2 


q(x,t) 


(B13) 


where 

S cross-sectional area 

T longitudinal tension in team 

El bending stiffness 

x longitudinal coordinate 

q transverse load 

w,v lateral and longitudinal displacements, respectively 

The term Ny/R occurs for the ring but not for the beam because the ring is 
curved . 


It should be noted at this point that the ring equations (eqs. (B9) and 
(BIO)) contain the same shortcomings inherent in the cylinder equations 
(eqs. (Bl) to (b 6)). For example, both sets of equations make the assumption 

that l/n.2 is negligible in comparison with unity (where n is the circum- 
ferential mode number) . As a result, these equations lose their accuracy for 
low values of n (e.g., n < 6). In the case of cylinders, Morley (ref. 25 ) 
removed this difficulty by modifying the bending terms. This modification was 

accomplished by replacing DV^w by D^V 2 + w in the final equations; the 

analogous modification of equations (B9) and (BIO) results in the following set 
of improved equations of motion for a ring (eqs. (la) and (lb)): 




St 2 


where 



+ 




+ 





+ ph ^ = q(y,t) 
St 2 


N-t; 




+ 



Except for the bending terms, equations (B9) and (BIO) and (la) and (lb) 
are identical. Because of this difference, however, only the latter set yields 
the correct linear ring vibration frequencies, and for this reason they are 
employed in the present analysis. 


31 


APPENDIX C 


APPLICATION OF THE METHOD OF AVERAGING 


The method of averaging (sometimes called "the method of slowly varying 
amplitude and phase") is often employed in nonlinear vibration problems. This 
technique is demonstrated here by applying it to equation (l^) which is repeated 
here for convenience: 



+ 2(3 c 





= G n cos Dt 


To apply the method of averaging to this equation, let 

£ c (t) = A(t)cos£dt + $(t)J 
= A cos X 


where A and $ are assumed to be slowly varying functions of t, and X des- 
ignates [h t + 0(t)J . Taking the derivative of £ (t) gives 

~ = -Afl sin X + — cos X - — A sin X (Cl) 

dT dT dT 


In the method of averaging, the result in equation (Cl) is replaced by two 
equations; namely, 

— - = -AD sin X (C2) 

dr 

and 

— cos X - — A sin X = 0 (Cj) 

dr d T 

The second derivative d 2 £ c /d.T 2 is then, computed from equation (C2), 

d 2i 

= _AD 2 cos X - — fl sin X - A — fl cos X (C4) 

dr 2 dT dT 

Next, the preceding approximations for £ c , d£ c /dT, and d^^dT^ are sub- 
stituted into equation (l^). The result is 
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(1 - n 2 ) 


Ik cos X - — Q sin X - A — ft cos X - 2j3 0 J2A sin X 
dr dT ° 


+ ^ A cos x(-G 2 A^cos 2 X - QA — sin X cos X - J2A 2 — cos^ 
2 V dT d T 


n 2 A 2 sin 2 X^ = Gn cos(X - $) 


(C5) 


Both sides of equation (C5) are then multiplied by cos X, and the result is 
added to equation (C3) after the latter has been multiplied by D, sin X. This 
procedure yields 

(l - H 2 )a eos^X - A 0, ^ - 2p c £lA sin X cos X + ^ A cos^ -fl 2 A 2 (cos^X - sin 2 x) 


-tokf— sin X cos X + A - cos^] 
\cLt dr / 


= G n (cos 2 X cos 0 + sin X cos X sin 3>) (C 6) 


Finally, this equation is "averaged" by integrating over one period on X. In 
the integration, A(t) and <&(t) are approximated by their average values A 
and $, for example. 



A(t)cos 2 X dX ~ 



A cos 2 X dX 


= Ait 




cos X dX 



= 52 A 5 
4 d T 


where d$/dT is the average value of d$/dT. When equation (c6) is averaged 
in this fashion, it becomes 

(-[ n 2)l onT efi 2 ^ 3eflA^ d$ n -r . 

VI - ft^/A - 2GA — r = G n cos $ (C 

dT 4 o dT 

In a similar fashion, a second equation for A and $ is obtained by 

(1) Multiplying both sides of equation (C5) by sin X 

(2) Adding this result to equation (Cj) after multiplying the latter by 
-0, cos X 

(3) Averaging the final equation by integrating both sides from X = 0 
to 2jr. 
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These manipulations give 


-20. — - 2j3 c flA - — = G n sin $ (c8) 

dT 8 dT 

Equations (C7) and (c8) are first-order coupled differential equations for 
A anrl $. They may be simplified to coupled algebraic equations when steady- 
state vibrations are considered. _For example, steady-state vibrations imply 
that the average values A and $ remain steady ( ij_e . , constant) with time. 

In this case, the average derivatives dk/dr and do/dT are identically zero, 
and equations (C7) and (c8) can be reduced to 

(l - o 2 )l - efl - A - = G n cos ¥ (l6a) 

-2p c HA = G n sin $ (l6b) 

These algebraic equations can now be solved simultaneously to give A and $ 
as functions of G n , e, p c , and 0. Such a result is indicated in 

equation ( 17 )* 
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STABILITY OF THE ONE-MODE RESPONSE 


To investigate the stability of the one-mode response (eqs. (l8)), both 
£ (t) and £ s (r) are disturbed slightly. The disturbed variables can be 

expressed as 

£ c (r) = A cos(flT + 3>) + ! c (t) (Dl) 


£ s (t) = 0 + 6 s (t) 


(D2) 


where I c (t) and | (r) represent small perturbations in the response of 

cos and sin respectively. The solution (eqs. (l8)) is said to be 

R R 

stable if £ c and £ s do not increase with time. 

Substituting equations (Dl) and (D2) for £ c and £ s into equations (12) 
and retaining only first-order terms in the perturbations gives 


and 


_2 

1 + cos 2 (f2r + $) 


* 2 *c eA 2 




J dT 


2 2 Sln 2(fiT + 


+ jl + |sin 2 (flr + $) - 2 cos 2 (Ht + $)jj £ ( 


dT 2 


2—2 

1 _ A COS 2(Qr + ¥) 


£c = 0 


(D3) 




where A cos(At + $) is presumed to satisfy equation (l^), and, for simplicity, 
the case p c = p s = 0 has been considered. (The case of nonzero’ damping is 
analyzed in ref . 15 . ) 


Analysis of Perturbation Equations 

It will be noted that equations (D3) and (DA ) are ordinary differential 
equations with coefficients that are periodic in time. Such equations arise 
frequently in stability analyses; they are typified by the Mathieu equation. 
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which has been studied extensively. (See refs. 21 and 26.) Its solutions are 
stable (bounded) for particular combinations of K and Q and unstable 
(unbounded) for others. Thus, associated with equation (D5) are certain stable 
and unstable regions in the KQ-plane. As indicated in the following para- 
graphs, the Mathieu equation results are useful in analyzing equations (D3) 
and (D^). 

As it stands, equation (Dj) has a periodic coefficient in every term. In 
simplified notation, it can be written as 


(l + p cos 2 x ) — ~ - 

d£„ 

pfi sin 2X — - + 

1 + |a(sin^X - 2 cos 2 x) 

V dr 2 

dr 

j 


(D6) 


where 



and X = ( + $ ) . 


The first-derivative term can be removed by using the following transfor- 
mation of variables: 


X = (At + <S) 

in = jL £LX 
dr dX dT 


= fl 


dX 




^X 


g c (X) = U(x)e 


psin2XdX 

2(l+pcos 2 x) 


(D7) 


= U(X)(l + p cos 2 x) 

With these new variables, equation (D 6) becomes 


J 


2 V ) d^U 
+ (JL cos^X) 


A 

ax 2 




+ p(sin 2 X - 2 cos 2 x) + p cos 2X + 


p 2 sin 2 2X - 
k (l + p cos 2 x) 


lU = 0 

(d8) 


Dividing by (l + p cos2x) and expanding equation (d 8) for small values of p 
gives the approximate Mathieu equation 



t-2, 


+ l) 


eA 2 


(fl 2 + l)cos 2X + 0 



U = 0 


where p has been replaced by its equivalent 


eA 2 


2~* 


(D9) 
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The first instability region of equation (D9) can be estimated from the 
known results for equation (D5). In terms of A and 0, the one-mode response 
is unstable (with respect to perturbations of within the area bounded by 


1 




< n < 1 


eA 2 



The boundaries of this region coincide with the locus of vertical tangents to 
the response curves, where the jump phenomenon of nonlinear vibrations occurs. 


In a similar fashion, the substitution 
tion (D4) to 


X = (fl T + $) 


transforms equa- 




cos 



0 


(DIO) 


The first instability region of this equation defines another_area in which the 
one-mode response is not stable. In terms of the amplitude A, instability of 
the solution requires that 


or 


j 2> 4(ai-x) 


for Q, 1 


a 2 > 


4(l - fl 2 ) 




for n 1 


These inequalities may be expressed as a restriction on H, and the solu- 
tion is unstable (with respect to perturbations in within the region 

1 - + o(e 2 ) < ft < 1 + + o(e 2 ) 

o o 

From a physical standpoint, the preceding conditions define a "critical 
amplitude" A cr of the single-mode response. To understand this result more 

fully, recall that ? c (t) is associated with cos ^ and that £ s (t) is 

related to sin 2^. The loading q(y, t) was chosen so that the driven mode 
R 

ny . ny 

cos receives excitation from the external force and the sin mode 
K R 

receives none. This difference in excitation shows up in the differential equa- 
tions which govern £ c and £ s , and a possible solution to these equations is 

£ c (t) = A cos(S2t + $) 

M t ) = 0 
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ny 

sin — mode does not vibrate. As pointed out by the stability 
R 

however, if the driven mode exceeds the critical amplitude 

A = H 1 - « 2 

Cr \ € S ! 2 

then the one-mode approximate solution becomes unstable. For amplitudes of 
cos —■ in excess of A cr , the companion mode will respond, despite the fact 
that it is not (directly) excited by the forcing function. 



where the 
analysis. 


Discussion of Stability Boundaries 

The first instability region and the response of the companion mode have a 

ny ny 

simple physical explanation. The cos — and sin — modes may be regarded 

R R 

as the only modes of vibration of a two-degree-of-freedom system, where both 
modes have the same natural frequency. As a result of the nonlinear terms in 
the equations of motion, dynamic coupling exists between the modes. When one 
ny 

mode cos is harmonically excited in the vicinity of its natural frequency, 

the dynamic coupling causes the companion mode to respond. (in the terminology 
of nonlinear vibrations, the companion mode is said to be "parametrically 
excited.") Many other vibration problems give rise to similar coupled-mode 
vibrations; one such example is the nonlinear sloshing of liquid in circular 
tanks . 

Up to this point, only the first (and most important) instability region 
has been discussed. The second, third, and higher instability regions of equa- 
tions (D9) and (DIO) indicate that the solution (eqs. (l8)) becomes unstable in 
narrow areas near 12 = l/2, l/j>, l/4, .... Additional calculations are 
required to provide adequate solutions in these regions; such results are dis- 
cussed on pages to 5 6 in reference 15 under the title "Ultra-Harmonic 
Response." (For a general discussion of ultraharmonic responses, stability of 
vibrations, and jump behavior, see ref. 26.) 

The results of the stability analysis for the one-mode response are sum- 
marized in figure 4, where the qualitative effect of damping is indicated by 
the dashed lines. Damping affects the stability boundaries primarily near 
Q = 1 and for small values of A. 
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CALCULATION OF A, B, $, AND t 


_ _Equations (22) and (23) form a set of five algebraic equations for A, B, 
<3>, A, and i|r. From these_equations, it is desirable to obtain response curves 
of the variations of A , B, <I>, and i)/ with 0 for constant values of 3 C , 

3s, e, and G n . These results are best obtained by using an indirect 
approach, as is frequently the case in nonlinear problems. 

The technique used was to treat G n as an unknown and compute B, G n , $, 
and i|r for several values of A at particular values of ft. This procedure 
resulted in curves of G n plotted against A with ft as a parameter. Then 

response curves of the variation of A with ft for constant values of G n 
were obtained_by cross plotting. In a similar fashion, curves of B plotted 
against ft, $ plotted against ft, and so forth, were determined. The method 
is outlined in the following section. Note that_in this_approach A and ft 
are regarded as being given, with B, G n , $, A, and i|r the unknowns . 


=52 


To compute G n , it was first necessary to represent sin 2A, cos 2A, and 


a in terms of known quantities by manipulating equations ( 23 c) and ( 23 d), 
For example, when B is nonzero, these equations become 


-2 

cos 2A = + 

A 2 


4 (ft 2 - l) 

eft^A 2 


(El) 


and 


_ 8 p a 
sin 2A = 


2-2 
eft A 


(E2) 


2 — 2 — 

respectively. By using the identity sin 2A + cos 2A = 1, these equations com- 
bine to give a quadratic for B 2 : 


^ , 8 (ft 2 - 1) -2 , 16 (n 2 - if + (8p s n) 2 _ (eft 2 ! 2 

_D + w ■** 


eft £ 


B + 


ur 


= 0 


*2 


Solving this equation yields two roots for B . In the present problem, the 
relevant expression is 


B 2 = 


_ 4(1 - ft 2 ) 

eft^ 


eft £ 


(eft 2 ! 2 ) - (8p s ft) 


1V2 


(E3) 
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The other root is similar in form but has a negative sign before the radical 
and was consequently discarded as an extraneous root because it will not satisfy 
equations ( 23 ) in the limit as 0 S tends to zero. 

— 2 

From this point on, B can be regarded as a known quantity, since the 
right-hand side of equation (E3) is expressed completely in terms of the given 

— O — 

variables. With B now completely determined, cos 2A can be found from 
equation (El). Similarly, sin 2A is given in terms of A, ft, e, and |3 S 
by equation (E2). 

From_an operational standpoint, _it is now a straightforward problem to 
compute B, A, G n , and <5. With A, ft, e, J3 S , and (3 C given, B 2 is 
found from equation (E3). Then sin £A and cos 2A are computed by using 
equations (El) and (E2). These results are all substituted into equations ( 23 a) 
and ( 23 b) which are repeated here for convenience 


G n cos 3> = (l - ft 2 ) A - ^^(a 2 - B 2 cos 2 a) 


G n sin $> = -20 c ftA - I^sin 2A 


and which relate the unknown quantities on the left to the known variables on 
the right. The nondimensional force G n and the phase $ can then be found 
from trigonometric identities; that is, 


G n = ^G n ^cos 2 <l> 


G^sin 2 # 


, 1/2 


(Ek) 


and 


$ = tan 


- 1 / u n 


G„ sin $ 


(G n cos $ 


(E5) 


where G n cos 0 and G n sin $ are computed from equations ( 23 a) and ( 23 b). 
Similarly, A can be determined as 


A = i tan^feS-^ 
\cos 2A/ 


(E 6 ) 


ince 
and 


sin 2A 
ijr are 


and cos 2A are 
easily obtained, 


known from equations (El) and (E2). 
since 


Finally, 
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and 

ijT = A + ® (E8) 

— 2 — — 

where B , A, and $> are all known from preceding calculations. 

This approach was used to obtain curves of the variations of A with G n , 

B with G n , and so forth for particular values of ft. By cross plotting the 
results, the response curves of figure 5 were obtained. The damping coeffi- 
cients, p s and 0 C , and the nonlinearity parameter e were held fixed during 

the calculations. The actual numerical work was done on an IBM 7090 electronic 
data processing system by using a simple FORTRAN program. 

By using a large value of 0 S in the calculations, the companion mode was 
suppressed and only the driven mode responded. This technique allowed the same 
program to_be used to compute both the one- and two-mode response curves. 

Plots of $ against ft and ijr against ft were also made, but they are not 
presented here since the phase of the response is of relatively minor interest. 
For the coupled-mode case, however, it_is interesting to note that A was 

usually close to zero. Both $ and i|r ranged from -n to -rt/2 as the fre- 

quency ratio varied from ft » 1 down to ft » 1, Similar results were obtained 
for the single-mode response. 


kl 


APPENDIX F 


STABILITY OF THE COUPLED-MODE RESPONSE 


As in the investigation of the stability of the single-mode response 
(appendix D), the stability of the coupled-mode response (eqs. (24)) may be 
examined by disturbing £ c and £ s directly: 

£ c (t) = A cos^Dt + $) + £ c (t) 

£ s (t) = B sin(flr + +) + | s (t) 

In the present case, however, this procedure results in two very complicated 
coupled differential equations involving £ c and | s , with each term in the 
equations having a periodic coefficient. A simpler analysis (which is adequate 
for purposes of this report) can be made by the method of averaging. 




Use of the Method of Averaging 

To apply the method of averaging to the stability problem, let 

A = A 0 + &(tj 
B = B 0 + b( t) 

$ = + cp(t) 

+ y(r) y 

where the steady~state solution has been redesignated by 

£ c ( t ) ~ ^-o ^os + ® 0 ) 

£ s (t) = B 0 sin(fi T + t 0 ) 


(FI) 


(F2) 


The variables a(-r), b(r), cp(t), and /(t) represent small perturbations in 
the amplitudes and phases of the steady-state response. The solution (eqs. (24) 
or (F2)) is said to be stable if a, b, cp, and 7 do not increase with time. 


Substituting the preceding expressions for A, B, $, and t (eqs. Fl)) 
into equations ( 21 ) and retaining only first-order terms in the perturbations 
gives coupled equations of the form 


da 

c-i .a + c 5 . — + c- .b + 
l<] dr 


C 4j 


db 

dT 


+ c 




+ c 


dcp , 

6 j S7 C 7J 


7 + c 


dT 


= 0 


(F3) 


42 



APPENDIX F 


The index j ranges from 1 to 4- and corresponds to each of the equations (21a) 
to (21d). The coefficients c^j are independent of t and depend upon such 

quantities as the original amplitudes, phases, and damping. For any set of 
steady-state values ^A OJ B 0 j H, e, etc.^ the values of cqj are constant, 

and equation (Fj) then represents a system of linear equations with constant 
coefficients . 

Such equations are readily solved by standard techniques. For example, 
when the substitutions 

a(-r) 
b( t) 

<P(t) 

7(t) 

are made in equation (F3), an eigenvalue problem for A results: 

|l - AM | = 0 (F4) 

where L and M are real, nonsymmetric four-by-four matrices involving the 
constants c ijj a^, bi, an< i 7q are a l so constants. When equation (f4) 

is multiplied out, it gives a polynomial which yields complex roots Aj_. If 

any root has a positive real_part, the perturbations increase with time and the 
steady-state solution (a 0 , B 0 > etc.) is unstable at that point. Con- 

versely, the vibrations are stable if all the roots Aq have negative real 
parts . 


= a 1 e 


.At 


= b n e 


At 


= 9i e 


.At 


7q e 


.At 


The Stability Boundaries and Their Interpretation 

For the case of no damping ^0 C =0, 0 S = o), equation (f 4) can readily be 
expanded to give 


+ 2a 2 (a 2 + 0 2 )a 2 + a 2 0 2 (a 2 - 0 2 ) = 0 

where the substitutions 

A = 2fiA 


(F5) 


p_ p 

2 _ eH Ap 


2 - 2 

O efi^B 0 
p 2 _ 
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have been used for simplification. The roots of equation (F5) will have nega- 
tive real parts if the inequality 


1 - 


2.38eA 0 2 

8 


+ 0 (e 2 ) < 2 


(F 6) 


is satisfied. 

On the other hand, in the absence of damping, equation (23c) shows that 
real nonzero values of B will exist only if 2 meets the condition 

2 ^ 1 + + o( e 2 ) (F7) 

o 

(Note that this condition agrees with that in the one-mode stability analysis.) 
(See appendix D.) The preceding inequalities combine to show that the approx- 
imate solution (eqs. (24)) is real valued and stable (in the absence of damping) 
within the region 

— P p 

1 _ g -^ £A + 0( e 2) < a < 1 + + 0 ( e 2 ) (F8) 

o 8 


The results of the coupled-mode stability analysis are presented in fig- 


ure 7* Real values of B do not exist to the right of the curve 2 = 1 + 


eA 2 

8 ’ 


and the two-mode solution is unstable to the left of the curve 2 = 1- 


2 . 38 eA 2 


Damping alters these stability boundaries, primarily near 2=1 and for small 
values of A. The effect of a small amount of damping is indicated qualita- 
tively in figure J by the dashed lines. 
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For reasons of conciseness and readability, it was necessary to omit fro: 
the body of the report many of the technical details pertaining to the experi 
ments. The purpose of this appendix is to explain some of these experimental 
intricacies which include the measurement of the forcing function and the ana 
ysis of the response by means of Lissajous figures. 

Equipment for Measuring Deflections 

Radial motions of the ring were measured by two inductance pickups. The 
deflection sensors were mounted on a fixture with a large bearing which allow 
them to move circumferentially around the ring. (See figs. 12 and 13 .) The 
signal from each pickup was fed through a carrier 
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amplifier and into one side of a differential amplifier, the output of which 
was in turn sent through a band-pass filter to a cathode-ray oscilloscope (CRO). 
(A block diagram of the arrangement is shown in fig. 8.) The system was con- 
nected so that the pickups operated in a push-pull fashion; this arrangement 
resulted in an antisymmetric operating characteristic, shown in figure 14- . 



Displacement, d, cm 

Figure 14.- Calibration curve for pickup system. 


Application and Measurement of Forcing Function 

Vibrations of the ring were excited by means of an electrodynamic shaker. 
The shaker was connected to the ring by a fine tungsten drive wire ( 0.001-inch 
(0.003-cm) diameter), as shown in figures 12 and 13 . The wire served as a very 
soft coupling spring between the shaker and the ring. Soft coupling was imper- 
ative, since the spring acted as a constraint and raised the effective natural 
frequency of the ring. 

The shaker amplitude and frequency were controlled by a standard 
oscillator-amplifier arrangement. By recording the displacement of the shaker 
and knowing the spring constant of the wire, it was possible to compute the 
force acting on the ring. That is, for a shaker displacement of S Q cos <nt, 
the amplitude of the force exerted on the ring was 

P = | k5 0 (Gl) 

where k is the spring constant of the drive wire, 6 0 is the amplitude of 
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the shaker displacement, and the factor l/2 arises from the geometry of the 
drive-wire arrangement. Strictly speaking, equation (Gl) applies only for 
static displacements of the coupling spring. However, when the fundamental 
resonant frequency of the coupling spring is much higher than the frequency 

supplied by the shaker, the equation P = i kS Q becomes a very good approxima- 
tion. In this case, the inertia of the coupling spring can be neglected, and 
the spring displacements are quasi-static. 

Although other types of coupling springs were tried, it was found that the 
tungsten drive wire was the only one which could satisfy the conflicting 
requirements of soft ooupling (i.e., a low value of k) and high resonant fre- 
quency. The wire used in the experiments had a fundamental frequency of 
300 cps (300 Hz) and a spring constant of only 0.035 lb/in. (6.12 N/m). 

Because of its low stiffness, the drive wire had only a slight influence 
on the natural frequencies of the ring. However, its attachment to the ring 

ny 

was such that the frequencies of the cos — modes were raised slightly, 

nv 

whereas those of the sin ~ modes were unaffected. This disparity created an 

K 

overlapping of two instability regions (fig. 15 ) and led to anomalous results. 
To counteract this effect, it was necessary to add a small concentrated mass to 

ny 

an antinode of cos which lowered its natural frequency and separated the 

R 

unstable regions. 

Adding the small mass had another beneficial effect, as it greatly reduced 

the tendency of the nodes of cos SL to "drift" circumferentially. Such 

R 

shifting of the nodes can be explained by considering small imperfections in 
the ring, as noted by 
Tobias (ref. 27). The 
added mass was a short 
piece of solder which was 
glued to the ring directly 
alongside the drive wire . 

This arrangement fixed the 
"preferential modes" of 
the ring. (See ref. 27.) 

The location of the 
mass and the drive wire 
may be thought of as the 
origin of the circumfer- 
ential coordinate y. In 
this case, the force 
exerted on the ring acts 
at y = 0. As a conse- 
quence, the cos —■ modes 
were driven by the forcing 


Y/'///A Driven mode unstable; "jumps ' 1 occur 
Pt-t-tSS Companion mode excited; both modes vibrate 




(a) l)q 5 (b) tj£ ijg. 

Figure 15.- Shift of instability regions due to changes in natural frequency. 
o c = Linear frequency of cos^ ; driven mode. u s = Linear frequency 
, . ny 

of sin companion mode. 
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function, but the sin ^ modes were not; that is, there was no generalized 

force on the sin modes. Note that the analysis was conducted for a similar 

K 

case. 

Actually, the loading used in the analysis (eq. (13)) was 

^ i 

q.(y,t) = cos x cos 

whereas the experimental loading was concentrated at y = 0 and can be repre- 
sented as 

q(y,t) = S(y)F cos cot (G2) 

In equation (G2), S(y) is the Dirac delta function and F is the load per 
unit length that the drive wire exerts on the ring. For the experimental 
arrangement described herein, F can be calculated from 

F = | (G3) 

where b is the length of the ring and P is given by equation (Gl). 

This type of loading (eq. (G2)) is analyzed in reference 15 . The results 
are very similar to those of the present analysis, provided that the driving 

frequency co is in the vicinity of the natural frequency of the cos ^ mode. 

In this case, the response of the ring can be approximated by using F in place 
of F n in equation ( 13 ) and other related equations. For example, the cor- 
responding nondimensional force Gn becomes 


G n = 


2 2 
itRph cdj^ 


(nRphb Jhcop^ 


(G4) 


where equation (G 3 ) has been used. 

A more detailed explanation of the experimental apparatus is given in ref- 
erence 15 ; as indicated therein, refined experimental techniques were used to 
minimize the nonlinearities introduced by the measuring system and by the sus- 
pension, shaker, and drive arrangement. 


Analysis of Response 

Detection of a single mode .- As noted previously, the pickup system sensed 
radial motions of the ring and converted them into electrical signals which 
were displayed on an oscilloscope. These signals were interpreted with the aid 
of Lissajous figures, as follows: 
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The voltage proportional to the radial displacement of the ring was fed 
into the vertical axis of the oscilloscope, while the horizontal axis was 
driven by the oscillator which controlled the shaker. This arrangement resulted 
in a Lissajous figure which indicated both the amplitude and the phase of the 
vibration that was being sensed by the pickup. By slowly moving the pickup 
around the circumference of the ring and noting the successive Lissajous fig- 
ures, it was possible to analyze the motion of the ring. For example, when the 

response consisted of primarily the cos — mode vibrating near its resonant 

frequency, successive 
Lissajous patterns such as 
those shown in figure 1 6 
were observed. The open ver- 
tical ellipses indicated that 
the major response was ±90° 
out of phase with the forcing 
function; measurements showed 
that the amplitude of the 

response varied as cos ^ 

around the circumference . 

At the zeros of cos — , the 

R 

Lissajous figures were small 
horizontal figure eights 
which indicated vibrations 
at twice the driving fre- 
quency. This pattern was in 
agreement with the analyti- 
cal results, which predicted the occurrence of this type of "double-frequency 
contraction . " 

These observations were in general agreement with the assumed vibration 
mode; namely. 


R 





Figure 16.- Circumferential variation of Lissajous figures. 
Single-mode response (resonance). 


w(y>t) = A n (t) cos ^ A n 2 (t) (G5) 

where A n (t) ~ cos cut. Subsequent measurements around the circumference of the 
ring supported this conclusion. 

Circumferential variation of response .- The voltage V from the pickup 
system can be represented as a Fourier series in time 

V(y,t) = V 1 (y)cos cut + Vg(y)cos Scot + V^(y)cos 3^>t + . . . (g6) 

with coefficients that vary around the circumference of the ring. The coef- 
ficients, V-^y) , V 2 (y)> . . ., contain the spatial variation of the deflection 

shape, combined with the harmonics introduced by the pickup system. By passing 
the signal through a harmonic analyzer and using a narrow 2-cps (2-Hz) band-pass 
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filter, the circumferential variation of the first, second, and third harmonics 
was determined. This determination was accomplished by exciting the n = 4 mode 
and recording the appropriate voltages at several positions along a half wave- 
length. The results are given in nondimen sional form in figure 17. The relative 
amplitudes of the harmonics were in the ratio 


As shown in figure 17, V-^(y) was almost a perfect sine wave. Similar 

results were obtained for other half waves of the n = 4 mode and led to the 

conclusion that Vj_(y) was very closely represented by mi r cos £*£■. It is 

believed that this voltage variation was due entirely to the motion of the ring, 
that is, 

w(y,t) ~ A cos cot cos ^ 

K 

would give rise to the observed results. 


The variation of the other harmonics was influenced by the nonlinearities 
of the pickup system which made it impossible to determine how much of the 
cos Scat term was caused by actual motion of the ring. As indicated in fig- 
ure 17, however, V^y) did contain a major component that was constant in 

space, as is suggested by the theory. 


The third harmonic in the signal resulted primarily from the nonlinearity 
of the pickup system. The data indicated that Vl(y) varied as 

Vj max sm-irtr, where r is the frac- 


tion of the half wavelength. Such 
a result can be explained by consid- 
ering the antisymmetric operating 
characteristic of the pickup system 
and noting that the major input to 
the system (i.e., the deflection of 
the ring) varied as sin jtr along a 
half wave. The latter result was 
confirmed by subsequent mode-shape 
measurements which showed that the 
primary deflection of the ring was 
independent of the vibration ampli- 
tude. (See fig. 9 and the related 
discussion in the body of this 
report . ) 

Additional results (ref. 15) 
include a plot of the root -mean- 
square response around the circum- 
ference of the ring and measurement 
of the amplitude of the "double- 
frequency contraction" at the nodes 



v. 


2.0 


1.0 £ 


2, node 


Second harmonic 


3, max 



Figure 17.- Spatial variation of response voltage. 
Amplitude of vibration: A max = 16.1h. 
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of cos — . The conclusion from all these measurements was that the experi- 
R 

mental deflection could be expressed as 

v(y,t) = A E (t)cos 22 + C E (t) (G7) 

when Lissajous figures of the type shown in figure 1 6 were observed. Amplitudes 

ny 

A E (t) and C E (t) are the experimental time-varying amplitudes of cos and 

the uniform contraction, respectively. The amplitude A E (t) was measured at the 

antinodes of cos 222, whereas the contraction Cg(t) was measured at the nodes. 

R 

Detection of two coupled modes . - Early in the experimental program, it was 
noticed that in some cases small changes in the driving frequency to caused 
significant modifications of the Lissajous patterns. This behavior was traced 

to the appearance of the companion mode sin 222 , which participated in the 

ny 

motion along with the cos r— mode. 


When coupled vibrations of this type occurred, successive Lissajous pat- 
terns such as those shown in figure l 8 were obtained. The open vertical 

ny 

ellipses were observed at the antinodes of cos ~ and indicated that the 

I\ 

response there was ±90° out of phase with the input force. At the nodes of 
cos 222, the Lissajous figures were straight lines, which is characteristic of 

a response directly in (or out of) phase with the input. The straight lines 
resulted from the sin 222 mode, which has its maximum amplitude at the nodes 

of cos 222. After trying various combinations of sines and cosines, it was 

found that the observed Lissajous figures correspond to a deflection of the 
form 

w(y,t) = A cos (cut + 3> A )cos 222 + b sin (cot + $ 3 ) sin 2^2 (g 8 ) 


which demonstrates the existence of the companion mode. 

When the amplitudes (A and B) and the phases ($ A and $ 3 ) of the motion 

are exactly equal, a pure traveling wave results. That is, equation (G 8 ) may 
then be written in the form 


■w(y,t) 


= A 


cos 


ny 

R 


(o>t + ® A ) 


(G9) 


Observations with the aid of a stroboscope disclosed the presence of traveling 
waves when Lissajous patterns such as those in figure 18 were noted. These 
observations gave further indication that the coupled vibration involved the 

cos mode and its companion mode sin 
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Figure 18.- Circumferential variation of Lissajous figures. Coupled-mode 
response (resonance). 

It was also noted that the appearance of the companion mode resulted in a 
reduction of the second harmonic content of the response. This result is in 
agreement with the analysis, as may be seen by examining the contraction term 

(eq. (9)) 

A 0 (t) = - £_[A n 2 (t) + B n 2 (t)] 

Since A n (t) and B n (t) are out of phase by 90° (i.e., A n ~ cos cot, 

B n ~ sin cut), the second harmonic component of A n (t) tends to cancel that 

p 

of B n (t) when the two terms are added. 

As a result of these and other observations, it was concluded that the 
experimental deflection could be expressed as 


w(y,t) = A E (t)cos + B E (t)sin ^ + C E (t) (G10) 

when Lissajous figures of the type shown in figure l8 were observed. The coef- 
ficients Ag(t), Bg(t) , and C-g,( t) are the experimental time-varying ampli- 
tudes of cos 22L, sin and the uniform contraction, respectively. 

R R 

Once the various responses of the ring had been identified, it was pos- 
sible to obtain response curves of the variations of A E (t) and B E (t) with 

co and to compare them with the analytical results . These tests are described 
in the body of the report . 
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APPENDIX G 


Other Re stilts of Interest 

In addition to the preceding results, several other experimental observa- 
tions were made. Part of this work was an experimental response survey, the 
results of which are given in table I. As indicated therein, several ultra- 
harmonic responses were noted, including some which involved coupling between 
two modes with different mode numbers. Results of this type are predicted by 
the infinite mode analysis of reference 15 . Additional experiments reported in 
reference 15 include the steady-state response of the n = 3 mode and some 
transient responses of the n = 4 mode. The latter responses demonstrate the 
growth of the companion mode with time when the driven mode exceeds the critical 
amplitude. Damping traces were also obtained for the n = 3 and n = 4 modes. 
For these and other details, see reference 15. 


53 



APPENDIX H 


ANALOG-COMPUTER STUDIES 


The approximate solutions obtained by the method of averaging indicate a 
gap in the coupled-mode -response curves. In order to find out what happens 
to the response in the gap region, an analog computer was used to simulate 
equations (12) which are repeated here for convenience. 


It 2 


2 Pc ~ +^c + |^c 


dT 




>C p 

dT^ 


c + I^C 


It 


+ i 


a 2 c 


S 2 

dr 


s + /45s 


dT 


= G n cos 


dT 2 


+ 2 Ps ^ + Cs + | 


>S dT 2 


&t 


+ 5 


d% 

C dT 2 


v dT 


= 0 


Analog- computer solutions were obtained for both the coupled-mode and 
single-mode responses. The procedure used and the results obtained are briefly 
discussed in the following paragraphs. 


In setting up the analog procedure, the parameters |3 C , (3 S , e, and G n 

were held constant, and responses were obtained at various values of Q. The 
procedure used was to set all the initial conditions equal to zero (i.e., 

\ 

£ c (0) = 0, S(0) - 0, etc .y and then turn on the forcing function at a par- 

ticular frequency ft. The driven mode would respond to the excitation, and 
usually it would reach a steady state, whether or not the companion mode 
also responded. Once steady- state conditions were attained, the stability of 
the response was examined by intentionally "bumping" or disturbing and £ s 

electrically. When the companion mode did not respond, disturbances to 


were found to damp out, as would be expected. Responses were obtained for 
two distinct cases which differ in the amount of damping in the companion 
mode . 


Case 1; coupled-mode response .- In the case of the coupled-mode response, 
the parameters used were 

Pc = Ps = °-002 (equal damping in both modes) 

e = 4.2 X 10~ 4 (nonlinearity) 

G n = 0.1 (force) 

These values are the same as those used to calculate the coupled-mode response 
by the method of averaging. As shown in figure 5, the analog results agreed 
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very well with the analytical solutions. In the gap region, the analog traces 
indicated nonsteady vibrations with beats. It was impossible to obtain steady- 
state responses within the gap, even though the problem was run for several 
hundred cycles. 

Case 2: single-mode response .- In the case of the single-mode response, 

the companion mode was heavily damped to prevent its vibration. The parameters 
used were the same as the preceding case except for p s which was p s = 0.05 

(companion mode damped heavily relative to the driven mode). The analog results 
verified the analytical solution for this case, as shown in figure 2. 
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TABLE I.- EXPERIMENTAL RESPONSE SURVEY 


n = 2 
n = 3 
n = 4 
n = 5 
n = 6 

n = 7 

n = 8 

n = 9 

n = 10 
n = 3 
n = 4 
n = 7 

n = 8 
n = 8 

n = 9 
n = 3 \ 
n = 5 / 
n = 4 l 

n = 6 ) 

n = 4 ) 
n = 7 / 

n = 4 \ 

n = 8 f 

n = 5 \ 
n = 9 ( 


Driving 
frequency, 
cps or Hz 

7-1 

17-5 

33-0 

53.2 

78.0 

107.4 

139.2 

179.6 

222.8 

8.8 

11.1 

213.0 

69.8 

47.6 

92.3 


39.0 


Type of response 




Resonance at driving 
frequency 


Ultraharmonic; order l/2 
Ultraharmonic; order l/3 
Subharmonic; order 2 
Ultraharmonic; order l/2 
Ultraharmonic; order l/3 
Ultraharmonic; order l/2 

| n = 3 at 26.7 cps (Hz) 

\ n = 5 at 53.4 cps (Hz) 

( n = 4 at 39.0 cps (Hz) 

\ n = 6 at 78 .O cps (Hz) 

1 n = 4 at 35.8 cps (Hz) 

( n = 7 at 107.4 cps (Hz) 

j n = 4 at 35.3 cps (Hz) 

( n = 8 at 141.2 cps (Hz) 
( n = 5 at 60.5 cps (Hz) 

\ n = 9 at 181.5 cps (Hz) 
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